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Abstract 

Expressions corresponding to the transmission of a uniaxial optically active crystal 
platelet are provided for an optical axis parallel and perpendicular to the plane 
of interface. The optical activity is taken into account by a consistent multipolar 
expansion of the crystal medium response due to the path of an electromagnetic 
wave. Numerical examples of the effect of the optical activity are given for quartz 
platelets of chosen thicknesses. The optical activity's effects on the variations of 
the transmission of quartz platelets as a function of the angle of incidence is also 
investigated. 
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1 Introduction 



Optical activity in crystals is a phenomenon which has been known for a long 
time (see Ref. [1] for an historical review). Two methods are currently used 
to describe the effect of the crystal optical activity on plane wave propaga- 
tion. The first one is phenomenological and is valid under normal incidence 
[2,3,4,5,6,7]. The second one is a first principle method and consists in solving 
the Maxwell equations with modified constitutive equations [1,8,9,10] {i.e. the 
relationship between the electric and magnetic fields vectors and the electric 
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displacement and magnetic induction field vectors). The choice of the consti- 
tutive equations is of prior importance [9] and was itself of phenomenological 
nature until the link between the optical activity as well as the magnetic- 
dipole and electric-octupole responses of the crystal medium was consistently 
formalised [11, 12]. This theory leads to constitutive equations different from 
those previously used. To the author's knowledge, the constitutive equations 
of Ref. [11] have not yet been used to calculate the transmission matrix of 
an optical active uniaxial platelet. They were only applied to describe the 
reflection of a plane wave by an interface between an isotropic achiral and a 
uniaxial chiral media [13], but for an optical axis parallel or perpendicular to 
the plane of incidence. It is the purpose of this article to provide this trans- 
mission matrix in a simple and usable form. We further restrict ourselves to 
two geometric configurations often encountered experimentally, namely the 
optical axis parallel and perpendicular to the plane of interface. 

Recent experimental results have shown that the contribution of the optical 
activity to the transmission matrix of a quarter wave plate can reach the per- 
cent level [6]. As an application of our transmission matrix formula, the origin 
of this effect is investigated and the case of half wave plates is also considered. 
Our formula being valid under oblique incidence, numerical calculations of 
the variations of the platelet transmission matrix as a function of the angle of 
incidence are presented. 

In section 2 we describe the method. In section 3, we derive the explicit expres- 
sions of the platelet transmission matrix for oblique and normal incidences. 
The contribution of the optical activity to the transmission matrix is studied 
numerically for various quartz platelet thicknesses in section 4. 



2 Formalism 

A monochromatic plane wave of wavelength A impinging a uniaxial optically 
active crystal slab of thickness I is considered. The crystal is assumed to be 
non-absorbing, non-magnetic and surrounded by an isotropic achiral ambient 
medium. 

A direct axis system x, y, z is defined such that the z axis is perpendicular 
to the plane of interface. The origin of the z axis (i.e. z = 0) is fixed on 
the first plane of interface so that z = £ corresponds to the second plane of 
interface (see Fig. 1). The unit vector basis attached to the system axes is 
denoted {x, y, z}. Without loss of generality, the plane of incidence is taken to 
be yz. To determine the transmission and reflection matrices of the interfaces 
between the ambient medium and the crystal faces, we follow the method of 
Ref. [14]. 
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Fig. 1. Schematic view of an optically active platelet with its optical axis in the 
plane of interface. The angles and axes used in the calculations are also drawn. 

Near the first interface z — 0, the electric and magnetic field vectors read: 



Incident: E = 
H = 

Reflected: E = 
H = 



A s s + A p p 
1 



exp(i[ujt — k ■ r]) 



-k x 



j2 UJ 

B s s + B p p 
1 



A s s + A p p 



exp(i[ut — k • r]) 



Refracted: E = 
H = 



/IqUJ 



(1) 

(2) 
(3) 

exp(i[cut - k' • r]) (4) 

exp(iut) (5) 
C (k 0+ x o + - in uTo + ) exp(-ik G+ • r)+ 



-k'x 



exp(i[ut — k' • r]) 
B s s + B p p' 



C Q o + exp(-ik Q+ • r) + C e e + exp(-ik e+ • r) 
1 



C e (k e+ x e + - ifi uTe + ) exp(-ik e+ • r) 



exp(iut) (6) 



and near the second interface z 



Reflected: E = 
H = 



C' Q o_ exp(— ik Q _ ■ r) + C' e e_ exp(— ik e _ ■ r) 
1 



[XqUJ 



exp(itot) (7) 
C' (k _ x o_ — ifj, uTo_) exp(— ik Q _ ■ r)+ 



Cg(k e _ x e_ — ifi LdTe_) exp(— ?k e _ ■ r) 



Refracted: E 
H 



A' s s + A' p p 
1 



exp(i[ut — k ■ r]) 



-k x 



A' s s + A'p 



exp(i[ut — k ■ r]) 



exp(iut) (8) 

(9) 
(10) 



with the incident field vectors given by Eqs. (5,6). 
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In these expressions, k and k' are the wave vectors in the ambient medium; the 
three vectors s, p and k form a direct basis, they are given by s = pxk/|pxk| 
and s = p' x k'/|p' x k'| with s perpendicular to the plane of incidence (i.e. 
s = x); T is a second-rank tensor coming from the constitutive equations [11], 
its components are related to those of the gyration tensor coefficients gij [2]; 
the electric vectors o± and e± are the solutions of the wave equation inside the 
crystal; the wave vectors inside the medium, k Q± and k e± , are given by the 
Snell's law and by the condition of existence of solutions of the wave equation. 

The wave equation inside the crystal is given by [11] 

k <* k P - S a(3 k 2 + fi uj 2 e a(3 + in u *VWy 1^ = (11) 

P \ 7 / 

where a,/3 and 7 stand for x or y or z and where the magnetic-dipole and 
electric-quadrupole response of the medium is embodied in the third-rank 
tensor A whose components are related to those of the gyration tensor[ll,13]. 
In Eq. (11), k a stands for the components of one of the four wave vectors k 0± 
and k e± ; Ep stands for one of the components of the four electric vectors o ± 
and e±; e a p are the elements of the dielectric tensor. 

The wave amplitudes are related by matrix relations. The refraction by the 
first (i.e. z = 0) and second (i.e. z = i) interfaces are described by [14]: 





t'pA °\ (12) 



V 

respectively. As for the reflections by the first and second crystal-ambient 
medium interfaces (inside the crystal), one gets 

C ^ R ( C '°) and ( C '°) = P:'RP + 1°°] (13) 



e / V e / V e/ V e ) 

respectively. In the previous equations, the following 2x2 phase matrices were 
introduced: 

(exp(— iip .) 
(14) 
exp(-^ e± )y 

with ip 0± = £(k 0± • z) and ip e± = £(k e ± • z). 



The 2x2 matrices T, T' and R are obtained from the continuity conditions 
of the projections of the electric and magnetic field vectors in the planes z = 
and z — £. Once known, these matrices are used to compute the platelet 
transmission matrix M such that E mit = ME in where E in and E mit are the 
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electric field vectors before and after the crystal slab respectively. Note that M, 
Ei n and E out are thus represented in the {§, p, k} basis, so that M reduces to 
a 2 x 2 matrix and E in and E out to two-dimension vectors when plane waves 
are considered. 

Taking into account the multiple reflections inside the crystal, one gets: 



M = T'P + T + T'P + RPZ l RP + T + ■■■ = T'P, 



1 - RP~ l RP+ 



T (15) 



where we have used Y^Lo X 1 = [1 — X] 1 and where 1 is the 2x2 identity 
matrix. 



We shall now choose a crystal symmetry class group in order to specify the 
two tensors A and T. 



3 Transmission matrix of a platelet with the optical axis parallel 
and perpendicular to the plane of interface 



Quartz crystal is widely used in crystallography to test the reliability of ex- 
perimental setups [7,6] and to manufacture retardation plates. We thus re- 
strict ourselves to the non-centrosymmetric crystals belonging to the symme- 
try classes 32, 422 and 622 [2,15] for which the tensors A and T have similar 
expressions [13]. In the crystallographic reference frame, i.e for the crystal 
optical axis aligned along the z axis, the non-vanishing components of these 
tensors are [13,15]: 

%x = % y = 933/2, f zz = g u - 5(33/2 (16) 

A X yZ AyxZ 933l Ay ZX A Z y X Qlll ^ZXy A XZ y Qll (^) 



with Ti m = fJ-ocTim and Ai mn = HocAi mn and where gu and 533 are the two 
independent coefficients of the gyration tensor [2]. Note that for another ori- 
entation of the optical axis, standard tensor transformations [2,15] are used 
to determine the representations of T and A. 



We shall now consider two geometrical configurations often encountered ex- 
perimentally: a platelet with its optical axis parallel and perpendicular to 
the plane of interface. The two cases are treated separately because different 
approximations are made during the calculations. 
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3. 1 Optical axis parallel to the plane of interface 



The wave equation (11) can be written in a matrix form 

(£ + 5£)E = (18) 

where the contribution of the optical activity is contained in the second term 
88 and where E = (E x , E y , E Z ) T . Writing the wave vector k = (u/c)(0, (3, k z ) T , 
with (3 = n a sin 9 and where 9 is the angle of incidence and n a the optical index 
of the ambient medium, we obtain: 



n 2 cos 2 + n 2 sin 2 — k z 
A sin cos 




A sin cos 
n 2 cos 2 + n 2 sin 2 — k 



7 2 

z 







and 





-iguh 



iguk z 





cos 2 + 5-33 sin 2 0) sin cos 0(533 -511) 



-i(3(gn cos 2 + 533 sin 2 0)^ 
i/3 sin cos 0(533 - 511) 







with A = n 2 — n 2 Q and where n Q = ^ e Q / eo and n e = \Je e /eo are the ordinary 
and extraordinary optical indices respectively. In the previous expressions, 
is the azimuth angle describing the orientation of the optical axis (note that 
= corresponds to the optical axis perpendicular to the plane of incidence). 



Since optical activity is induced by the electric-quadrupole and magnetic- 
dipole responses of the crystal medium[ll], we assume that \gij\ <^ \n e — n \. 
The condition of existence of a solution for Eq. (18) reads det(£ + 5S) = 0. 
This condition gives us the four possible values for k z : k oz± = ±(&£^ + Sk oz ) 
and k ez± = ±(&^ + Sk ez ) which correspond to the ordinary and extraordinary 
waves respectively (the sign ± refers to the direction of propagation). The 
zero order terms kf) and kf) are given by the solution of det(£) = 0: 



n\ -/f (cos 2 0+ ^fsin 2 0) 



1/2 



(19) 
(20) 



6 



and the first non-vanishing contributions in g^ read as 



Sk n 



5k P 



with 



r,2 



Anlk 



(9n n o - I 3 sin ^(#33 ~ 9u) ~ Sk) 



(fl'ii^o - ft sin2 4>{9zz ~ 9iif + Sk) 



(21) 
(22) 



5k 



7,2 



[nt 



13 sin"0)A 



4 2/2, 2\ 

n 9n{n e + n ) 



+(3 n sin 4>[2g n (g 33 - g n )(n 



n. 



rA 



+(3 sin 4 0(^3 -giiYK + ni) 



As expected [11], Eqs. (21) and (22) are of second order in g^. 



(23) 



Accordingly, the ordinary and extraordinary electric field vectors are decom- 
posed as follows: o ± = o±^ + 5o± and e± = + 5e±. The zero order terms 
and are solutions of £o±^ = and Se^ = 0: 



(o) 
M 



N (hWsm<i>,-kW cos<P,±P 



COS i 



ef = N e {n 2 cos 0, {k^) 2 sin 0, T jfeW sin 0) T . 



(24) 
(25) 



1 and |e±^| 



where N Q and N e are the normalisation factors such that |o± 
1. 5o± and 5e±, the first order contributions in gij/A, are the solutions of 
£5o± = —8£o±* and E5e± = —SSe^. Although k oz ± and k ez ± are of second 
order in g^, implying that det(S) = at first order, it turns out that the 
conditions of existence of solutions for these equations are fulfilled. We obtain 



5o± 



±iN n 



5e 



± 



~ 2 

cos <f)[gun 2 + $ sin 2 (p(g 33 - g u )] 
sin (f)[g u n 2 + $ sin 2 <p(g 33 - g^)] 
■P sin <p[gun 2 + g 33 A + ft sin 2 <f)(g 33 - g u )]/(n 2 - ft) 

1 ~kf) sin <Pn 2 /(n 2 - ft) 
—kf) cos0 



~2 

±iN e [g 11 n 2 Q + (3 sin 2 (f>(g 33 - g u )] 



/?cos0p(° 2 )) 2 -A]/(n 2 -/5' 



(26) 
\ 

(27) 



Once the ordinary and extraordinary wave and electric vectors are known, the 
matrices 71, T and T are determined by supplying the boundary conditions 



7 



in the planes z = and z — £. For an oblique incidence, these calculations 
are more simply performed numerically. Compact analytical expressions are 
obtained for the special case of normal incidence and will be given bellow. 



To define a platelet transmission matrix at first order in g^'/A, we introduce 
the three matrices R , T and Tq obtained by setting g n = g 33 = to zero in 
the expressions of R, T and T respectively. We then define the matrices SR, 
and ST and ST' such that R = R + 5R, T = T + ST and T = + ST. 
Having thus isolated the zero order from the first order terms in g^/A, we 
define the platelet transmission matrix M = M + SM and thereby 



Eout = E<2t + 5E out such E^ t = M E in , SB out = 5MB in . (28) 



,(o) 



Proceeding as in Eq. (15), we obtain: 



SM =ST'{T^)- 1 M + M Tq 1 ST+ 



T' 



T' 



1 — P + R P_ 1 R 
1 — P+RqP- 1 Rq 



P + RqPZ 1 SR 



n -l 



1 — P+RqP- 1 Rq 



P+T + 



P+SR 



1 — P_ 1 RqP + Rq 



PZ l RoP + T G (29) 



and 



1 — RqP- 1 RqP+ 



-i 



(30) 



Note that M only depends on the optical activity through the phase matrices 
P + and Pi 1 . 



Under normal incidence, one gets k oz± = ±n [l — gf l /(2A)}, k ez± = ±n e [l + 
gf 1 /(2A)] to the first order in g u /y/A and the corresponding electric vectors 



i ■ a _l m o9n cos r 
°± — ( sm ± : , — COS (J) ± 



e± = (cos ± 



A 

in e gw sin . 



rnpgn sm (p T 

A ' / ' 



, sin =p 



tn e gu cos < 



-,of. 



(31) 
(32) 
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The reflection and transmission matrices are: 



R 



T 



in, 



—in 



n +n a 
2giiK-rag)+ff 3 3A 
° (na+n )(n a +n e )A 

2gu(n e n a +nl)+g 33 A 



2g 11 (n 2 -n 2 a )+g 33 A ' 
'■ (n a +n )(n a +n e )A 



n e +ria 



a (n a +no)(n a +n e )A 

2Wg 



n e +n a 



m a 



T' = 7Z((j)) 



in 



2gn (n e n a +n' 2 )+g 33 A 
" (n a +n )(n a +n e )A 

-2n 



~2n a 
rio+ria 
2grLno{n a +no)+gy i A 
(n a +n )(n a +n e )A 

2n e 



rie+ria 



n +n a 



-in 



2giino(na+no)+g 33 A 
e (n a +n )(n a +n e )A J 



(33) 
(34) 
(35) 



where 7Z((p) is the 2x2 matrix representing a rotation of an angle around 
the z axis. 



In these expressions, terms in g n are enhanced, unlike those in g 33 , by a factor 
l/(n e — n Q ) ^> 1. This is a major difference with the case of an optical axis 
perpendicular to the plane of interface (see next section). From Eqs. (26-27) 
one also sees that the contribution of #33 to the transmission matrix becomes 
noticeable only under oblique incidence. 



Let us mention that Eq. (34) agrees with the result of Ref. [9] for an isotropic 
chirality, i.e. gn = #33. It is however in disagreement with the expression of 
Ref. [10] obtained using phenomenological constitutive equations. The expres- 
sion of our platelet transmission matrix therefore disagrees with the one of 
Refs. [17,18]. 



The zero and first order platelet transmission matrices are given by Eqs. (30) 
and (29) respectively. For Eq. (30), we obtain the following simple expression[19]: 



/ 



M = n{(j>) 



4n n a exp(—iip ) 



(n a +n ) 2 -(n a -n ) 2 exp(-2iip D ) 





4re e n a exp(— iip e ) 



(n a +n e ) 2 -(n a -n e ) 2 exp(-2iip e ) 




with <f = (u/c)lk oz+ and ip e = (u/c)£k, 
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3.2 Optical axis perpendicular to the plane of interface 



When the optical axis is parallel to the z axis, one can write Eq. (11) as 
follows: 



K - P igaakz -iguP 



V 



E„ 



0. 



\ E 'J 



-ig^k z n 2 -k; (3k z 

ignP Ph n\ - (3 2 J 

To first order in g^, the four solutions for k z are given by k oz± — ±(o — b) 1 / 2 
and k ez± = ±(a + b) 1 / 2 with 

/3 2 



a = n °-2n 2 



X + ni) 



(37) 



b = Ugh + %rWi - 3^33) + J£[2(n* + n 2 e )(g 33 - 9ll ) 2 + (n 



1/2 



,2\2i 



(38) 



and the corresponding electric vectors read as [16] 

/ 

v = N 



n 2 e (nl-k 2 )-n 2 J 2 " 



ik z g^n 2 e + i(3 2 k z (g 33 - g n ) 
-0[k 2 z (g33 - 9u) + 9unl] 



where v = o±, e± for k z = k oz± , k ez± respectively and where N is the normal- 
isation factor such that |v| = 1. 

As in the previous case, the transmission matrix is more simply computed 
numerically. Under normal incidence, one gets k oz± = ±(n Q —g 33 /2) and k ez± = 
±(™o + #33/2) to the first order in g 33 . The corresponding electric vectors are 
the two circular polarisation states that we write = (=pi, 1,0)/ \/2 and 



3 ± 



= (±i, 1,0) /\/2. For the transmission and reflection matrices, one gets: 



V2n a I i 1 



,T> 




n n - n n 



1 



n a + n. 



n Q + n a \ iq 



Using Eq. (15), the following platelet transmission matrix is obtained 

/ 



M 



2n n a [exp(-iy? ) + exp(-iy? e )] 
(n a + n D ) 2 - (n a - n Q ) 2 exp(-i[<p + ip e \) 



1 -tan((</? e -<p )/2) 

\tan(((p e -(p )/2) 1 

(39) 
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with <p = £k oz+ and ip e = £k ez+ . This expression agrees with the result ob- 
tained using a 4 x 4 matrix method [20]. 

Note that the interface matrices of eq. (39) do not depend explicitly of the 
coefficients gu and #33. Thereby, the variations of M with these coefficients 
come only from the phases tp and tp e . Note as well that the coefficient gu 
contributes only under oblique incidence. 

4 Numerical results 

To estimate the effects of the gyration coefficients on the platelet transmission, 
an experience closed to the HAUP[7] is considered: a linearly polarised He-Ne 
laser beam crosses a quartz platelet and then a supposed perfect analyser. The 
incident electric vector reads E in = (cosip, sin-0) T in the {§, p} basis (i.e. in 
the basis {x, y} under normal incidence, see Fig. 1). After the analyser, whose 
eigen axis corresponds to, say, the x axis, the beam intensity is measured as 
a function of ip and 0. 

For our wavelength, A = 0.6328/iin, the optical indices are [21] n = 1.542637 
and n e = 1.551646 and the gyration coefficients g n = —5.9 • 10~ 5 and g 33 = 
10.1 • 10~ 5 [7]. 

To first order in g^/A, the beam intensity after the analyser is given by: 

I x = |E<W • x| 2 + [(E£ ■ x)(5E out ■ *)• + (E2 ■ x)*(5E out ■ *)] (40) 
where the field vectors are defined in Eq. (28). 

We first consider the normal incidence and a crystal optical axis parallel to 
the plane of interface. As for the plate thickness, two values are chosen: l\ = 
87.80/iin and £ 2 = 1000.93/iin (i.e. a first and a fourteen order quarter wave 
plates). 

The intensity I x is shown as a function of ip and in Fig. 2 (it is similar for 
the two plate thicknesses i\ and £ 2 ). Defining I 0x = |x ■ E^l 2 = |x • M E in | 2 , 
the difference A x = I x — I 0x is shown as a function of and i\) for the plate 
thickness i\ in Fig. 3. The effect due to the medium gyrotropy is surprisingly 
high for ip ps 7r/4 and the same behaviour is indeed observed for the second 
plate thickness i 2 - In the expression of Iq x , the contribution of the optical 
activity comes only from the phases f and ip e . 

To distinguish the chiral contribution coming from the phases ip Q and ip e from 
the one coming from the interface matrices, we define Iq x = |x- MoEi n | 2 where 
M is given by Eq. 36, but fixing ip = 2nn /\ and ip e = 2nn e /\ (i.e. by 
neglecting the optical activity). The difference A 0x = I 0x — I 0x is shown in 
Fig. 4 as a function of i\) and for the plate thicknesses l\. These figures 
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indicate that the large variations observed in Fig. 3 come from the chiral 
dependence of the interface matrices. For the plate thickness £2, the shape of 
A 0x is identical to Fig. 4 but scaled by a factor of ten. As expected, the chiral 
dependence induced by the phase shift alone becomes only sizeable for thick 
plates. 




Fig. 2. I x , for a normal incidence angle and for the plate thickness £±, as a function of 
the angles (j) and ip (in radian). The optical axis is parallel to the plane of interface. 




Fig. 3. A x , for a normal incidence angle and for the plate thickness £\, as a function 
of the angles <fi and ip (in radian). The optical axis is parallel to the plane of interface. 




Fig. 4. Aoa;, for a normal incidence angle and for the plate thickness £±, as a function 
of the angles <j) and ip (in radian). The optical axis is parallel to the plane of interface. 



The previous discussion holds for quarter wave plates. Let us consider another 
plate thickness £3 = 1018.49/xm, i.e. a fourteen order half wave plate. I x , A x 
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and Aox are shown as a function of and ip in Figs. 5, 6 and 7 respectively. 
The effect induced by the optical activity is here reduced by two orders of 
magnitude with respect to a quarter wave plate and is dominated by the 
variation of the phases ip and ip e with g n . For a first order half wave plate 
(£4 = 105.36/xm), the effect on A 0x is reduced by an order of magnitude 
whereas A x is essentially unchanged. 




Fig. 5. I x , for a normal incidence angle and for the plate thickness £3, as a function 
of the angles 4> and ip(in radian). The optical axis is parallel to the plane of interface. 




Fig. 6. A x , for a normal incidence angle and for the plate thickness ^3, as a function 
of the angles 4> and ip (in radian). The optical axis is parallel to the plane of interface. 




-5e-05 
-0.0001 
-0.00015 



Fig. 7. Aoa , for a normal incidence angle and for the plate thickness £3 , as a function 
of the angles (f) and ip (in radian). The optical axis is parallel to the plane of interface. 
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As mentioned in section 3.1, under normal incidence and for the optical axis 
parallel to the plane of interface, the contribution of the optical activity to 
the platelet transmission matrix comes essentially from the coefficient gn. The 
contribution of the coefficient #33 appears under oblique incidence and can, a 
priori, be distinguished from the other one by considering various angles of 
incidence. We now consider the ratio r = I s /I p where I s and I p are obtained 
by substituting s and p to x in Eq. (40) respectively. To reproduce the ex- 
perimental results of Ref. [6], we also fix ip = 7r/4, i.e. E in = (1, 1) t /a/2. The 
ratio r is shown as a function of (ft and 9 for the plate thickness l\ in Fig. 
8. To distinguish the contributions from gn and (733, we show in Figs. 9 and 
10 the relative differences 5r = (r — tq)/t and S'r = (r — r{)/r where ro and 
T\ are obtained by setting gu = g 33 = and gn = in the expression of 
r respectively. From these figures, one sees that, as reported in Ref. [6], the 
contribution of the optical activity reach the percent level for certain values 

of (f). 




Fig. 8. r as a function of the angles 9 and 4> (in radian) for an incident electric 
vector Ei n = (1> ' /V% and for the plate thickness t\. The optical axis is parallel 
to the plane of interface. 




Fig. 9. 5r as a function of the angles 8 and 4> (in radian) for an incident electric 
vector Ej n = (1> 1) T an d for the plate thickness i\. The optical axis is parallel 
to the plane of interface. 
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Fig. 10. S'r as a function of the angles 9 and 4> (in radian) for an incident electric 
vector Ei n = (1, l) T /\/2 and for the plate thickness l\. The optical axis is parallel 
to the plane of interface. 

In addition, we observe here that the contribution of #33 is negligible under 
normal incidence but reaches the few 10 -3 level for reasonable values of the 
angle of incidence. The shape of the contribution of gu and (733 being different, 
it seems possible to determine both coefficients by measuring r as a function 
of (f) for various angles of incidence. 

We finally consider the case of an optical axis perpendicular to the plane of 
interface. We choose £5 = 1000/im for the plate thickness and, as previously, 
Ein = (l)l) T /v / 2 for the incident electric vector. In addition to Sr, we in- 
troduce the phase angle a defined by the argument of the complex number 
(E ou t • s)/ (E ou t • p) and the difference Aa = a — «o where «o is obtained by 
setting gu = in the expression of E out . Figs. 11 and 12 show 5r and Aa as 
a function of the angle of incidence 8. 
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-0.016 
-0.018 

0.05 0.1 0.15 0.2 0.25 
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Fig. 11. 5r as a function of the angle 9 (in radian) for an incident electric vector 
Ej n = (lj I) 7 '/v2 and for the plate thickness £5. The optical axis is perpendicular 
to the plane of interface. 
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Fig. 12. Act, in radians, as a function of the angle 8 (in radian)for an incident 
electric vector Ei n = (1, l) T /\/2 and for the plate thickness £5. The optical axis is 
perpendicular to the plane of interface. 

As in the case of an optical axis parallel to the plane of interface, the mea- 
surements of r and Act as function of the angle of incidence lead to the de- 
termination of both (?n and (733. The sensitivity to these parameters is even 
higher when the optical axis is perpendicular to the plane of interface. 



5 Conclusion 

The transmission matrix for a uniaxial and chiral platelet has been calcu- 
lated for an optical axis parallel and perpendicular to the plane of interface. 
The optical activity has been taken into account by considering the electric- 
quadrupole and magnetic-dipole responses of the crystal medium. Simple and 
usable expressions have been derived for normal incidence and for the optical 
axis parallel and perpendicular to the plane of interface. The case of an oblique 
incidence has been treated partly analytically and partly numerically. 

As an application, numerical studies of the transmission of quartz platelets of 
various thicknesses have been provided. Under normal incidence, it was shown 
that the contribution of the optical activity can reach the percent level for 
quarter wave plates, independently of the plate's order. Such a high contribu- 
tion was found to be due to the transmission and reflection interface matrices. 
The case of half wave plates was also investigated and, unlike quarter wave 
plates, the main contribution of the optical activity came from the phase shift 
induced by the optical path inside the crystal. The variation of the intensity 
after a quartz platelet was studied as a function of the angle of incidence and 
of the orientation of the optical axis. A sensitivity to both gyration coeffi- 
cients gn and g 33 was observed under oblique incidence. This suggests that 
varying the angle of incidence one could determine experimentally these two 
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coefficients with a unique crystal cut. With respect to the results of Ref. [22], 
oblique incidence also offers the possibility to determine very accurately the 
platelet thickness and thereby the gyration coefficients. 
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